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A study of eccentricity (en) fiuctuations and its possible impact on final state momen¬ 
tum anisotropy for symmetric collisions are presented in the framework of Glauber 
model. Effect of fiuctuations of nucleon’s position on the initial geometry has been stud¬ 
ied using a new method, where the difference between oppositely moving spectators is 
taken as a measurement of eccentricity fluctuations. This study shows that higher har¬ 
monics (n =3, 4 and 5) of eccentricity are less sensitive to fluctuations in transverse plane 
compared to the 2^^ harmonic. Position fluctuations in transverse plane will increase £2 
and hence possibly V 2 for the most central nucleus-nucleus collisions. For semi-central 
and peripheral collisions, the fluctuations have opposite effect, it deceases the eccentric¬ 
ity £ 2 - The fluctuation of initial geometry can be studied in collider experiments by 
studying the spectator distribution on the both sides of the beam. 

PACS-key : 25.75.Ld, 25.75.-q 


1. Introduction 

One of the main goals of the high energy heavy-ion collision experiments is to 
study the QCD phase diagram i>2,3,4^ achieve this goal, one has to understand 
the properties of the system formed in such collisions. The momentum azimuthal 
angular anisotropy parameter has been considered as a good tool for studying 
the system formed in the early stages of high energy collisions at Super Proton Syn¬ 
chrotron (SPS), Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider 
(LHC) 5,6,7,8,9,10,11,12^ describes the harmonic coefficient of the azimuthal 
Fourier decomposition of the momentum distribution with respect to the reaction 
plane angle (d>) The final state momentum anisotropy (vn) reflects the hydro- 
dynamic response of initial spatial anisotropy (e„). According to hydrodynamical 
description, Vn is sensitive to the geometry of initial state of the system formed 
in the collision as well as the hydrodynamic evolution governed by the equation of 
state of the matter 8,9,10,11,12,14^ 

Knowing the initial geometry and fluctuations in heavy-ion collisions has recently 
been shown to have important consequences on interpreting the experimental data 
from various experiment at RHIC and LHC. Experimentally measured non zero 
odd harmonic (n>3) has been interpreted as the result of statistical fluctuations in 
the transverse positions (according to uncertainty principle) of nucleons undergo¬ 
ing hadronic scattering. Moreover, measured V2 cannot be described by an smooth 
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initial energy density distribution, unless one includes flow fluctuations arising due 
to the eccentricity fluctuations in the calculations addition to Vn, there 

are other experimental observables which cannot be explained without including 
eccentricity fluctuations. For example, dihadron correlations in azimuthal angle 
and pseudorapidity The contribution from the odd harmonics associated with 
the particle azimuthal angle distribution to dihadron correlations is found to be an 
important factor. 

Several phenomenological studies have been carried out on initial geometry anisotropy 
and fluctuations to understand its influence on experimental data 20,21,22,23,24,25^ 
The aim of this paper is to discuss the centrality dependence of various harmonics 
of initial spatial anisotropy and its sensitivity to the fluctuation in position of nu¬ 
cleons using a new method. In this paper, £„ are calculated within a framework 
of Monte Carlo Glauber (MCG) model which allows the generation of colli¬ 
sions with event-by-event-fluctuating initial condition. The paper is organized in 
the following way. In section 2, Glauber model has been briefly discussed. Section 
3 describes the study of e„ and its sensitivity to the fluctuations using the MGG 
model. Finally, the summary has been given in section 4. 

2. Model Description 

In MGG model, the nuclear distribution function inside a nucleus is taken to be 
of the Woods-Saxon form, 

^ l + exp[{r-R)/dy 

where the radius (i?) and the diffuse constant (d) are taken as R = 6.38 fm, d = 
0.535 fm for Au nucleus. In this model, nuclei are assembled by positioning the 
nucleons randomly in a three-dimensional coordinate system, on an event-by-event 
basis, according to the Woods-Saxon density profile. A collision between two nuclei 
is considered as a sequence of independent nucleon-nucleon collisions. In a nucleus- 
nucleus collision, two nucleons with transverse distance d < yJcfNN/ t^ will collide 
with each other. Here aMN is the total nucleon-nucleon cross-section. 

Using the transverse position coordinates of each colliding nucleon, various mo¬ 
ments of participant eccentricity have been calculated as: 

-y/(r2 COs{ruppart))‘^ + SUY{rHPpart))‘^ 

where r and ippart are the polar coordinate positions of participating nucleons. 

r = \/x^ + y'^ (3) 

and 

(Ppart = tan~^{y/x). (4) 

In this study, approximately 7 million events for each configuration with fixed 
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Figure 1: (Color online) Distribution of nucleons in the transverse plane for one 
typical Au+Au event at 5 = 5 fm. Solid and open circle represents the participants 
and spectators nucleons, respectively, from two colliding nuclei. 

impact parameter 6 = 1 /m, 5 fm, and 9 fm, are generated for Au+Au collisions 
with unn = 40 mb. Variation in (Tnn does not change results qualitatively. A 
standard MCG model code which is used as an input in AMPT model has 
been used to generate events in this study. The distribution of nucleons in the 
transverse plane for a single event is shown in Fig. 1. 

3. Results and Discussion 

Fig. 2 shows magnitude of spatial initial eccentricity in transverse plane for 
different harmonics (from n=2 to n=5) in Au+Au collisions at 6 = 1 fm, 6 = 5 fm 
and 6 = 9 fm . The large value of £2 for 6=5 fm and 6 = 9 fm is due to initial 
elliptic shape of the overlapping region in a collision of large impact parameter. 
All odd higher harmonics (n > 2) are generated due to fluctuations in transverse 
positions of nucleons. For a nucleus with smooth density distribution, all odd higher 
harmonics (n > 2) of will be zero. In case of nucleus-nucleus collisions at 6 = 
1 fm, the initial overlapping geometry is almost isotropic, hence magnitude of £2 
is small and comparable with values of £3, £4 and £5. Centrality dependence of 
£2 can be understood, since it reflects the anisotropy of overlapping region of two 
nuclei. But we observed, as shown in Fig. 2, that all higher harmonics of £„ reveal 
similar centrality dependence, like to that £2. The fluctuations behave differently for 
collisions with different impact parameter. Although the non-zero £2 is originated 
due to initial elliptic shape, it can be modified by the nucleon density fluctuations. 
Therefore, it is very important to understand the role of initial geometry fluctuations 
in heavy-ion collisions. The main purpose of this paper is to investigate how various 
harmonics of eccentricity change with nucleon density fluctuations. 
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Figure 2: (Color online) Magnitude of spatial eccentricity in transverse plane for 
different harmonics in Au+Au collisions at 6 = 1 fm, 6 = 5 fm and 6 = 9 fm. 


The fluctuation varies on event-by-event basis. Therefore, the number of spectators 
(S) and participants (Npart) will also vary on event-by-event basis. Moreover, in a 
single event, number of spectators from target nucleus (labeled as A) and projectile 
nucleus (labeled as B) can be different due to the fluctuation of nucleon’s position 
in the transverse plane In this study, the difference in the number of spectator 
between two colliding nuclei (|5'^ — S’®!) has been used to quantify the fluctuation, 
more fluctuation means large difference and vice-versa. Total number of events for 
a hxed impact parameter are divided in several sub-groups based on | — 5''® |. 

Fig. 3 shows, event-by-event distribution of number of spectators in A-nucleus {S"^) 
and in B-nucleus(S'®) from MCG model at 6 = 5 fm. The maximum difference 
between and can be of the order of 50. 

Values of £ 2 , £ 3 , £ 4 , and £5 as function of \S^ — S'-®] are shown in Fig. 4. Panel 
(a), (b) and (c) corresponds to events with fixed 6 =1, 5 and 9 fm, respectively. 
In each case, magnitude £ 3 , £ 4 , and £5 are scaled to match their values with £2 at 
Pqj. central events with 6 =1 fm (i.e. panel(a)), all harmonics of 
£„ increases with increase in \S^ — S'^|. This indicates that, fluctuations enhance 
the anisotropy for the most central collisions. As we have fixed the value of impact 
parameter, changes in £„ is entirely due to fluctuations. What is striking in this 
observation is that, fluctuations making the system more elliptic, while also makes 
the system more triangular, quadratic and pentagonal. 

Now for semi-central events with 6 =5 fm (panel(b)), we can see that there is 
small change in eccentricity for n > 3, their values are increasing with increase in 
.^g observed that the magnitude of £2 decreases with increase in 
unlike central events. This shows that the fluctuations in the transverse 
plane are decreasing initial elliptical geometry a nucleus-nucleus collision. 

For peripheral events with 6 = 9 fm (panel(c)), £2 changes sharply with change 
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Figure 3: (Color online) Number of spectators in A-nucleus {S^) vs Number of 
spectators in B-nucleus {S^) in Au+Au collisions at 6 = 5 fm. 


in — S^\, and shows a decreasing trend with increasing fluctuations, like semi¬ 
central events. The 63 and £5 increases with fluctuations and almost negligible 
change for £4 with respect to other harmonics. 

To quantify the sensitivity of £„ to fluctuations, ratios between maximum 
and minimum (£™“) value of eccentricity has been calculated. The maximum 
change in £ 2 , £ 3 , £4 and £5 are ~ 15%, 12%, 3% & 7.5%; ^ 8.3%, 4.5%, 2% & 3.5%; 
and 52%, 14%, 3% & 8.5%; for b =1, 5 and 9 fm, respectively. This indicates 
that of all the harmonics, £2 is more sensitive to the fluctuation and then £ 3 . This 
observation is consistent with the previous study done using different asymmetric 
collision in AMPT model where it was shown that V 2 is more sensitive than 
V3. On the other hand, we observed that the £4 is less sensitive to the fluctuations 
compared to £ 5 . 

We know that the final state momentum anisotropies are driven by the initial spatial 
anisotropy and flow coefficients (u„) are proportional to £„. Therefore the change 
in En due to fluctuation will affect in similar manner. Sensitivity of the 
to the fluctuation can be different compared to £„ and that depends how the £„ 
evolve through different stages of the fireball history and translate into final-particle 
momentum anisotropies. But qualitatively, sensitivity of the Vn could be similar 
like £„. Therefore, from Fig. 4 we expect that the initial fluctuations in transverse 
plane will generate more V 2 in most central collisions, whereas for semi-central to 
peripheral collision magnitude of V 2 will be reduced due to the fluctuations. This 
observation is in agreement with previous study carried out using 3-1-1 Viscous 
Hydrodynamics in Ref 

In real collider experiment, one can easily measure number of spectator in both 
direction of beam and hence the difference between them. Therefore, using this 
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Figure 4: (Color online) Spatial eccentricity £ 2 , £ 3 , £ 4 , and £5 versus \S^ — S''®| for 
collisions at (a) & = 1 fm, (b) 6 = 5 fm and (c) b = 9 fm 


method one can identify the events with different amount of fluctuations in one 
centrality bin and can measure Vn to understand the effect of fluctuation. Narrow 
centrality bins will be more appropriate for this study. In experiment, centrality 
is usually estimated by number of particles produced in the collision. The number 
of spectators is anti-correlated with the number of particle participating nucleons 
in the collisions. We have studied how the average Npart changes with 
as shown in Fig. 5, to estimate the effect of event mixing from different centrality 
group. We can see from Fig. 5 that the change in < Npart > is less than 3%. 


4. Summary 

A study on initial collision geometry fluctuations for a symmetric system us¬ 
ing MCG model has been presented. It has been observed that all other higher 
harmonics of £„ show centrality dependence like £ 2 . A new method using number 
of spectator nucleons has been used to separate events with different amount of 
fluctuations. Due to fluctuations in the transverse plane of colliding nuclei, £„ (and 
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Figure 5: (Color online) Average Npart versus |5'"^ — S'^| for collisions at 5 = 1 / to , 
b = 5 fm and 6 = 9 /to. Here < Npart > for collisions at 6 = 1 /to and 5 /to are 
scaled by factor 0.28 and 0.4, respectively. 

possibly v„) increases in most central collision. For semi-central and peripheral col¬ 
lisions, 62 (and possibly ^ 2 ) is minimised by the fluctuations, on the other hand all 
other higher harmonics are found to be higher due the fluctuations in the transverse 
plane. Moreover, we observed that 2"^^ harmonic is more sensitive to the collision 
geometry fluctuation compared to higher harmonics although higher harmonics are 
generated due to fluctuations. This new proposed method can be applied to more 
realistic transport model and in real experiment to study the fluctuations in 
which is very crucial to understand various properties like transport coefficient of 
the system created in heavy-ion collisions. 

Only the fluctuation of nucleon’s position in the transverse plane has been discussed 
in this paper. For ultra-relativistic nucleus-nucleus collisions, nuclei are contracted 
along beam axis (Z-axis) and looks like as thin plates in lab-frame. As a result, 
fluctuations of nucleons position along longitudinal directions are negligibly small. 
But in case of collisions at a low energy, like AGS energy, where the out-of-plane 
squeeze-out phenomena in elliptic flow was observed, fluctuations along longitudi¬ 
nal directions may not be negligibly small. Future investigation can be done in this 
direction using various transport model to see the effect of longitudinal fluctuations 
in the final states momentum anisotropies using this new proposed methods. 
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